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SUMMARY 


A numerical analysis is presented for the problem of a flat rectangular 
rubber membrane with a circular rigid inclusion undergoing high strains due 
to the action of an axial load. The neo-hookean constitutive equations are in- 
troduced into the general purpose TITUS program by means of equivalent nookean 
constants and initial strains. The convergence is achieved after a few itera- 
tions. IhB method is not limited to any specific program. The results are in 
good agreement with those of a Company sponsored photoelastic stress analysis 
The theoretical and experimental deformed shapes also agree very closely with 
one another. For high strains it is dsmonstrated that using the conventional 
HOOKE law the stress concentration factor obtained is unreliable in the case 
of rubberlike material. 


INTRODUCTION 


The structure of a radial motor vehicle tire is made up of two types of 
components namely the reinforcing cords and the rubber. The most immediate 
problem It. tire 3tress analysis is that of the large displacements in the in- 
flated tire described in a previous paper (Reference 1). It appears that the 
mo**s important components are the reinforcing cor^s allowing the tire to take 
a stable inflated shape. This particular problem can now be considered as 
solved. 

However, in order to solve the complete problem, the rubber's behavior 
must also be adequatly analyzed by means of an as economical as possible modi- 
fication of existing programs. Up until now, this vary challenging problem of 
non linear material behavior and incompressibility has only been solved in a 
few special cases (Reference 2). 

The aim of the present work is to stress analyze the rubber parts of the tire 
by using NASTRAN and TITUS. A test specimen encompassing a rigid inclusion is 
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culculuLnd In uniaxial uxtnnslnn and the results 
phnlnulunLi c Lty . The experimental evidence shows 
these methods now available to the designer. 


are compared with those nl 
tin: limits of validity of 


SYMBOLS 


( K i Stiffness matrix, Nm 

(0! Matrix relating strain to nodal displacements, mm/mm 

A Surface of the membrane element, m 2 

CD) Constitutive law in matrix form, Nm ' 2 

Stress, Nm" 2 
£ Strain, mm/mrr, 

S 0 Initial strain, mm/mm 

V POISSON'S ratio, (no units) 

W Clastic potential per unit volume of the unstrained body, Nm" 2 

and C 2 Constants of MOONEY, Nm' 2 

Ii Strain invariants, i = 1,2 and 3, (no units) 

(E) Neo-hookean constitutive law of a membrane in matrix form, Nm" 2 

Initial hydrostatic stress of a rubber membrane, Nm' 2 
u Displacement, m 

(Tj and 0*2 Principal stresses in the middle plane of the membrane, Nm" 2 

^ Radius of curvature of the transverse isostatic, m 

s Curvilinear abscissa, m 

01 Normal stress tangent to the edge of the disc, Nm" 2 

C Photonlastic material constant, N'^m^ 

a Radius of the disc, m 

Subscripts: 

T transposed 

t true 

y coordinate perpendicular to the loaJ axi 1 ; centered i.n the middle of 

the inclusion 

y coordinate along the load axis centered in the middle of the inclu- 

sion 



GENERAL APPROACH 


HOOKE'S LAW 


The elementary tests carried out show that for the material under consi- 
deration HOGKE's law applies in relating true stress * to strain even in the 
3b mm/mm range. Thus, it may be assumed tentatively that the non linearity 
of the rubber's constitutive law is only the result of the large displacements 
experienced by the rubber. 


Stiffness matrix 

The stiffness matrix of the membrane plate element used can be written 
(Reference 3) : 

(k) = ( B 1 j x (0) x (B) x A x d 

d Thickness, in 

A Surface of the membrane element, m 2 

(B) Matrix relating strain to displacement, mm/mm 

(0) Constitutive law, Nm~2 

As a result of the incompressibility condition Ad = const. The accuracy 
of the forces and the displacements depends on the accuracy of the terms R and 
D. 

Definition of 0 



Experimentally the uniaxial law is 

(Tj s Eg where E is YOUNG'S modulus 

As a result <T must be replaced by in the equations. The only remain- 
ing term that has to be calculated in the course of the extension of the spe- 
cimen is B. Thus, this problem would seem to be identical to that of the large 
displacement problem. 


* The true stress is computed per unit section area of the deformed body 
whereas the conventional stress is computed per unit section area of the 
undeformed body. 


51 


In the case of rubber, however, it is well known that there is an addi- 
tional. pressure type, term "p" in the constitutive equation. Tn eliminate 
“p". represented by (PI jl,) above, an additional equation in terms nf dis- 

placements must be written for each element. NASTRAN with its scalar element, 
can handle such an equation. The resulting data input problem is however very 
cumbersome. Thus, this solution may be uneconomical for every day use. 

In order to demonstrate the existence of "p" a large displacement calcu- 
lation was carried out with £. = 0 and v> = .5. The largest transverse dis- 

placement along the x-axis passing through the middle of the inclusion was 
13 % in error relative to the experimental values. This error was much larger 
than the one obtained by means of the theory developed below. 


MOONEY-RIVLIN CONSTITUTIVE LAW 


The most common type of rubber material behavior equation is that of 
flooney-Rivlin (Reference 4). Considering that W is the elastic potential mea 
sured per unit volume of the unstrained body the postulated function is 

W = C 1 CI 1 - 3) ♦ C 2 (I 2 - 3) 

1^ are the strain invariants (i = 1,2) 

C., C n are the constants postulated by Mooney. 

I C 


The theory of plane stress of very thin membranes applies to the rubber 
specimen considered here. The deformations are symmetric about the middle 
plane of the body and are essentially uniform throughout the thickness. The 
pressure type component ”p" is eliminated because in the present problem the 
normal stress perpendicular to the specimen's surface is zero. Large displace- 
ment equations are used in the b matrix. 


The equations obtained are: (O’) * (E) {£.} ♦ 


(T are defined at points in the deformed body, but are measured per unit 
area of the undeformed body. 

The E and C T 0 are functions of not only and C£ but also of £ X x* Eyy ar >d £ X y* 
The matrix E is positive definite in the strain range considered. 


In uniaxial extension the above equations in terms of true stresses must be 
identical to the well known equation (Reference 5) s 



and 



1 + e„ 



1 

1-v 2 


where V = ,5 


This happens only if C2 is zero and the de‘ ,-mations are limited in size. Such 
a material is called neo-hookean. The constant C-j is determined by means of 
the latter equations in an uniaxial elementary extension test. C 1 * . 71 MN m" 2 . 
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The solution of the equations is carried out by using an equivalent se- 
cant modulus method. The full load is applied in this trial and error approach. 
The first solution in obtained by the hookean constitutive low where <J 0 is set 
at /mu. 


NASTRAN AND TITUS ANALYSIS 

The program of J.T.Qden, OKI, (Reference 6), is designed for analyzing 
rubberlike material. Thus it came under consideration first, but it can only 
analyze plane strain plates whereas our problem is a plane stress problem * . 

The solutions obtained by the TITUS and NASTRAN programs have been compa- 
red at the first iteration. TITUS uses isoparametric quadrilateral membrane 
elements while NASTRAN has constant stress CQDMEM elements. The stresses differ 
only by 2 However the difference between the displacements of t .aSTFIAN as 
compared to those of TITUS was 4 %. The results of NASTRAN were much further 
away from the experimental ones than those of TITUS. In thin particular case 
the CPU computation time was SO s for TITUS and 84 s for NASTRAN using UNIVAC 
1108 (EXEC 8). 


MODIFICATION OF TITUS 


The TITUS program was developed in France by CITRA now called SPIEBATI- 
GNOLLES Inc. Because of the pi'oximity of the development team it was easy to 
modify the program. By means of a minor modification it is possible to calcu- 
late the modulus E and <r o internally elementwise at each iteration with the 
help of C-j and the strains. 

The test of convergence was carried out by comparing the arithmetic mean 
of the displacements obtained at each iteration. At lower leads (9.8 - 19.6 N) 
the convergence was achieved after about six iterations whereas at 29.4 N ten 
iterations were needed. In the first case the computation time was 84 s CPU 
on UNIVAC 1108 (Exec 8). 


MODELING OF THE PLATE 


The finite element idealization of the membrane encompasses 107 nodal 
points and 84 quadrilateral elements. In order that the theoretical solution 
and experimental results could be satisfactorily compared, the three loading 


* In linear elasticity plane strain and plane stress problems are conjugate. 
This is not the case, however, for rubberlike materials. 
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conditions were given in terms uf displacements at the end of the specimen. 
Ihe modeling of the membrane is shown in Fig. 1. 


Boundary conditions 
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Sines the loading and the deformations are assumed to be symmetrical, only 
one-quarter of the plate needed to be considered. 


EXPERIMENTAL WORK 
TEST SPECIMEN 


The model test specimen is a rectangular coupon cut out of a polyurethane 
plate furnished by PHOTOLASTIC Inc. The coupon is tran- parent and isotropic 
when not loaded. A circular hole is cut out of its cent^ ' and is fillet, ’n 
with araldlte which is reinforced with glass beads at a ratio of 100 ' . Tne 
stresses due to the contraction of the disc during polymerisation have bean 
observed by means of crossed polarisors and have been eliminated by an applied 
compression load in order to keep the neutral state of stress in the specime: . 
The disc is much sti x fer than the rest of the coupon and there is perfect ad- 
hesion between them. The grips are glued on to the ends of th8 rectangle. The 
only load applied is a vertical load along the specimen's axis and it is mea- 
sured by means of strain gages. Viscoelastic effects are suppressed by load- 
ing up gradually. 


The dimensions of the specimen are 117 mm x 42 
meter of the disc is 14 mm. 


x 1.02 mm and the dia- 
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Using a transmission polariscope, the isoclines are determined with the 
help of in-plane polarized white light and the isochromes with the hslp of mo 
nochromatic circularly polarized light (0 = 5690 A). 

The isostatics are obtained from the isoclines using graphical means. 


PHOTOELASTIC STRESS ANALYSIS 


In rubber the lightwave path difference is proportional to the difference 
of the principal stresses CT-) - 0"2 and also proportional to the instantaneous 
thickness of the specimen (Reference 7) 

$ * c e (cr, - <t 2 ) 

The material constant C is determined by a uniaxial elementary tension test. 
The value obtained is 

C * 3.21 ♦ 0.03 m 2 daN' 1 


The principal stresses along the vertical and horizontal symmetry axes are de- 
termined by integrating graphically the equation of Lam6 and Maxwell (Reference 


d, -Oi 


-t- 

e, 

— <■ m 0 

2 s 1 

a, - Oi 

*Sl . 0 
ts* 


where f> is the radius of curvature of the transverse lsostatic and s Is the 
curvilinear abscissa at a given point. The subscripts 1 and 2 ref er tc the 
two families of isostatics. 

The starting point of the integration along the x-axis is taken at the 
edge of the specimen where the stress <T2 is zero. 

For the integration along the vertical axis the poin of reference for the 
integration is taken in the region of uniform stress betwe. n the grips and the 
disc where d 2 is zero. Along the edge of the disc the stresses (T-j and 0*2 are 
obtained using the normal stress tangent to the disc and they vary as fol- 

lows: 

* o; _ "6 CT rt 

i ~T- ■ Tr- 

end CT, » O; -*■ (<T, - <0 Sin 2 <x 

where a is the angle between the direction of Q and the isostatic CT -j . 


The value of <T, is a function of the accuracy of the measurement of the 
Isoclines. As the experimental determination of the latter is relatively inaccu- 
rate. in particular at the top of the disc, the accumulated errors may be quite 
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large. Tor thin reason, starting from the horizontal axis, 
and cr ^ at the top aie 1 5 % greater than the ones obtained 
other axis. 


the values of cr-) 
starting form the 


DISCUSSION OF THE EXPERIMENTAL AND THEORETICAL RESULTS 

DEFORMED SHAPE 

The theoretical and experimental results obtained for the deformed shape 
are in excellent agreement as shown (Figs. 2-4). By different experimental me- 
thods the overall transverse displacements at the horizontal symmetry line have 
been determined as follows: 


LOADING 

EXPERIMENT 

TITUS 


u x 

u x 


(mm) 

(mm) 

6 * 

- 0.4 ♦ 0.05 

- 0.434 

11.5 % 

- 0.75 ♦ 0.05 

- 0.806 

20.5 % 

- 1.4 ♦ 0.05 

- 1.44 


ISOSTATICS 


The theoretical and experimental results showing the distribution of the 
isostatics over the surface of the rubber coupon are plotted in Fig. 5. On the 
left side are shown the calculated principal stresses and on the right side 
the envelopes of the corresponding experimental principal stresses. Taking into 
account the fact that the theoretical results ore relative to the undeformed 
surface of the specimen, the agreement is again excellent. The following table 
shows the values of the applied longitudinal force. 


L0A0ING 

EXPERIMENT 

TITUS 


Force 

Force 


(N) 

(N) 

6 H 

9.8 ♦ 0.2 

10.6 

11.5 % 

19.6 ♦ 0.4 

18.4 

20.5 \ 

29.4 ♦ 0.6 

28.2 


The mesh used in modeling the ends of the specimen was very coarse, the prin- 
cipal aim being to demonstrate the behavior of an inclusion imbedded in a rubbe 
matrix. Thus, the error obtained is accordingly larger in this region. 
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ISOCHROMES 


The Hxpurlmnntiil strum; distribution over Lhu surl.icu ul' Lliu sptH.lmuns 
has huun ustablishmJ liy moans of the isochromes. No quant J tot. Lvn compnrlsnn I:. 
shown horn on the numerical results are relative to tho untie Formed shape and 
are in terms of the conventional stresses, (Fig. 6-8). The automatic plotting 
of the required true stress isochromes is being developed at the present time. 

However, the shape of the isochromes shown agrees qualitatively with the 
experimental ones. As demonstrated below the numerically obtained longitudinal 
true stress concentration factor is very accurate. Thus it can be conjectured 
that the agreement must be also quantitative. 


PRINCIPAL MEMBRANE STRESSES 


In Fig. 9 the true principal membrane stresses together with the experi- 
mental ones are shown. The shape of the two families of curves obtained are 
identical. However there is a vertical shift of the theoretical ones relative 
to the experimental ones. The difference is quite small and remains within the 
limits uf the accuracy of the experiment. It must be noted that at the top of 
the inclusion, on the y-axis, the experimental results differ according to 
whether the point is approached from the right or the left. The mean of the 
two values is located very near to the theoretical point. 

At the intersection of the x-axis with the contour of the inclusion two 
nearly identical compression stresses are obtained experimentally. This re- 
sult agrees with those obtained by theoretical consideration in reference 9. 
The numerical calculation gives two stresses of opposite sign, however. This 
is explained by the fact that the stresses are calculated at the center of 
gravity of the element. In this region C of Fig. 10, the stress gradient is 
very large. Thus, even though the 0T^ stress is positive at j = 1.07, the cal- 
culation point, it is negative at • 1. that is at the experimental recording 
point. Taking these facts into account the agreement between the finite ele- 
ment results and those of the experiment is very good. 


STRESS CONCENTRATION FACTOR 


I he stress concentration of the longitudinal principal stress along the 
y-axis is plotted in Fig. 10. The maximum stress concentration factors are 

Experiment 1.28 

Finite element results 

Linear classical elasticity <1*54 


The agreement between the experiment and the numerical results is excel- 
lent. It can be concluded further that the llnrar elasticity gives unreliable 
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values in the case of rubber for strains reaching 1U \ ui more. 


CONCLUSION 


The close agreement obtained between theoretical and experimental results 
demonstrates the validity of the large displacement equations, and of the neo- 
hookean constitutive law used in the modified TITUS program. However, the use 
of the derived method is not limited to any specific program. After some minor 
modifications any geometrically nonlinear finite element program may bn ap- 
plied to the analysis of rubber at relatively high strain. 

The importance of using the proposed theory instead of the conventional 
HOOKE typa formulation to design rubber parts is made evident by the fact 
that, using the conventional theory, the stress concentration factor of the in- 
clusion obtained has an error of about 20 
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TITUS Solution 
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Experiment 


Figure 3. * Displacement Vectors, 


Load ■ 9.81 Newtons. 
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Experiment 


Figure *. - Displacement Vectora. Load » 29.4 Newtona. 
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Conventional Principal Stresses 
on the Undeformed Shape 

(TITUS) 


Experimental Isostatxcs 
- Principal True Stress Envelopes 
on the Deformed Shape 


Figure 5. - Conventional Principal Stresses and Experimental Isostatics 
Load ■ 19.6 Newtons. 


64 



TITUS Solution 


Experiment 


Difference of the 
Conventional Principal 
Membrane Stresses 
in the Undeformed Body 


Isochromes - Difference of the 
True Principal 
Membrane Stresses 
in the Deformed Body 


Figure 6. - Isochrom6s - Numerical Results and Experiment. 
Load ■ 3.61 Newtons. 
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TITUS Solution 

Difference of the 
Conventional Principal 
Memorane Stresses 
in the Undeformed Body 


Experiment 

Isochromes - Difference of the 
True Principal 
Membrane Stresses 
in the Deformed Body 


Figure 7. - Isochromes - Numerical Results and Experiment . 
Load ■ 19.6 Newtons. 


TITUS Solution 


Difference of tha 
Conventional Principal 
Membrane Stresses 
in the Undeformed Body 


Experiment 


Isochromos ■ Difference of the 
True Principal 
Membrane Stresses 
in the Deformed Body 


Figure 8. 


Iaochromes - Numerical Results and Experiment 
Load » 29.4 Newtons. 









